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Abstract. We study the instability of magnetic helds in a neutron star core driven by the 
parity violating part of the electron-nucleon interaction in the Standard Model. Assuming a 
seed field of the order 10^^ G, that is a common value for pulsars, one obtains its amplification 
due to such a novel mechanism by about five orders of magnitude, up to 10^^ G, at time scales 
~ (10^ — 10^) yr. This effect is suggested to be a possible explanation of the origin of the 
strongest magnetic fields observed in magnetars. The growth of a seed magnetic field energy 
density is stipulated by the corresponding growth of the magnetic helicity density due to the 
presence of the anomalous electric current in the Maxwell equation. Such an anomaly is the 
sum of the two competitive effects: (i) the chiral magnetic effect driven by the difference of 
chemical potentials for the right and left handed massless electrons and (ii) constant chiral 
electroweak electron-nucleon interaction term, which has the polarization origin and depends 
on the constant neutron density in a neutron star core. The remarkable issue for the decisive 
role of the magnetic helicity evolution in the suggested mechanism is the arbitrariness of 
an initial magnetic helicity including the case of non-helical fields from the beginning. The 
tendency of the magnetic helicity density to the maximal helicity case at large evolution times 
provides the growth of a seed magnetic field to the strongest magnetic fields in astrophysics. 

Keywords: neutron stars, cosmic magnetic helds theory, magnetic helds 


Contents 


1 Introduction 1 

2 Magnetic helicity evolution in a cooling NS core 3 

2.1 Initial conditions 6 

2.2 Growth of electric conductivity for a cooling non-superfluid NS 6 

2.3 Magnetic helicity growth driven by weak eN interaction 8 

2.4 Change of the sign for CME 9 

3 Generation of magnetic fields in magnetars driven by weak eN interactions 10 

4 Discussion 11 

A Averaged electron current induced by chiral effects 13 


1 Introduction 

The study of strong magnetic fields inherent in some compact astrophysical objects, like neu¬ 
tron stars, is one of the hottest topics of the modern astrophysics [1]. This research received 
a significant impact in the wake of the observations of anomalous X-ray pulsars (AXP) [2] 
and soft gamma-ray repeaters (SGR) [3]. Some observational characteristics of SGRs and 
AXPs distinguish them from the more common accretion-powered pulsars in massive X-ray 
binaries and place them into a separate class of astrophysical objects. 

From the point of view of modern astrophysics [4], SGRs and AXPs are supposed to 
be highly magnetized, B > 10^® G, neutron stars, or magnetars. It should be noted that a 
typical pulsar can have a magnetic field up to Bq 10^^ G. The magnetic field Bq can be 
generated in the core of a supernova (SN) progenitor at some earlier stage of its evolution, 
and subsequently amplified during core collapse, for example, simply by the flux conservation 
[5]. The source of the strong magnetic field in a magnetar, in its turn, is still unclear. 

It was shown in refs. [6, 7] that the most popular model of the magnetic field in a 
magnetar, based on the turbulent dynamo [8], confronts with some of the observational data, 
and thus should be corrected. Some other models based, e.g., on the idea of the fossil field, 
are also discussed in the literature [7]. 

Anyway, if one adopts a dynamo scenario for the creation of the magnetic field in a 
magnetar, one has to propose a mechanism to amplify the seed field Bq = 10^^ G, which is 
typical for a pulsar, by at least three orders of magnitude. Recently, in ref. [9], we put forward 
a new mechanism which can provide this kind of amplification. That scenario is based on the 
magnetic field instability because of the parity violation in the Standard Model (SM) for the 
electroweak electron-nucleon (eN) interaction of ultrarelativistic degenerate electrons with 
non-relativistic degenerate nucleons in a neutron star (NS) as the SN remnant and a future 
magnetar. 

On the first glance the mechanism proposed in ref. [9] resembles the results of ref. [10], 
where the chiral magnetic effect (GME) in an external magnetic field, which is proportional to 
the chemical potential of massless charged fermions, was established. The contribution of the 
electroweak interaction into the averaged electron current was discarded in the case of massive 
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fermions. That result of ref. [10] is based on the ambiguous renormalization procedure for the 
photon polarization operator calculated in the first order of the perturbation theory in Gp, 
where Gp ~ 1.17 x 10“® GeV~^ is the Fermi constant. Nevertheless, the possibility that a 
non-zero weak interaction term can appear in other models has not been ruled out in ref. [10]. 

To avoid nonsequential method in ref. [10] during calculations of the electric current 
driven by electroweak interactions, in ref. [9] we considered the additive eN interaction 
in the same Dirac equation for chiral states of massless electrons in an external uniform 
magnetic held B. We obtained that the total induced electric current J, which is additive to 
the standard ohmic current Johm in plasma, entering the Maxwell equation, 

J = ^(1^5 + ^5)B, (1.1) 

TT 

is given by the difference of the right and left electrons chemical potentials, //s = (/iR —^l)/ 2, 
and includes the non-zero weak interaction coefficient 

U = ^ _ (1 _ , (1.2) 

where nn,p are the densities of neutrons and protons inside NS, ^ = sin^ = 0.23 is the 
Weinberg parameter and a^m ~ (137)“^ is the hne structure constant in eq. (1.1). 

We would like to stress the important role of the magnetic helicity density h{t) = H/V = 

J (A • B)d^x in the generation of magnetic helds B = V x A in magnetars^. 

Despite the problem of the magnetic held evolution in a SN remnant was treated in a 
simplihed way in ref. [9], we could derive some of the characteristics of the magnetic held in 
a magnetar, which are close to the observed ones. Nevertheless we should mention some of 
the shortcomings of our model in ref. [9]. Firstly, a stationary electric conductivity cJcond = 
const ~ (Tg)”^ at the typical temperature T ~ 10® K for a cooling neutron star instead 
of its decreasing temperature Tg(t) = r(t)/10®K was used in the analysis. Secondly, we 
assumed the maximal and monochromatic helicity density spectrum: h{k,t) = 2p-Q{k,t)/k^ 
and h{k,t) = h{t)6{k — k^), where pB{k,t) is the spectrum of the magnetic energy density. 
Thirdly, we considered the magnetic held with the largest scale Ar = k^^ = i^NS = 10 km 
only. 

Of course, the canonical inequality h{k,t) < 2pB{k,t)/k (see, e.g., ref. [11]) with a 
running wave number k, kmin < k < fc ma v. for the continuous helicity density spectrum h{k, t) 
and the magnetic energy density spectrum pB{k, t), as well as the diminishing temperature for 
a cooling NS core, dT/dt < 0, are more realistic conditions for our MagnetoHydroDynamical 
(MHD) problem, and we can improve our model in ref. [9] using all these ingredients. 

Our work is organized as follows. In the main section 2 we consider the magnetic helicity 
evolution in a cooling NS core starting from the Faraday equation generalized in SM due to 
the anomalous electric current in eq. (1.1). Basing on the Faraday equation and using the 
Adler anomaly, we derive in section 2 the full set of self-consistent kinetic equations for the 
three functions: the magnetic helicity density spectrum h{k,t), the magnetic energy density 
spectrum pB{k,t) and the CME parameter p^{t). After formulation of the initial conditions 

^It is well-known that, in a perfectly conducting fluid with fixed boundary conditions, the magnetic helicity 
H is a, conserved quantity. Even for a finite electric conductivity the magnetic helicity changes slowly. It leads 
to the existence of stable equilibrium configurations for magnetic fields. To the extent that the helicity is 
conserved, a non-equilibrium or unstable magnetic held with a finite conductivity cannot decay completely, 
since the helicity of a vanishing field is zero. 
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in section 2 . 1 , we explain in section 2.2 our choice of the initial time for which the law of 
the temperature cooling within a NS core should be valid. Then we present our results with 
the numerical solutions of the three kinetic equations for the magnetic helicity density h{t) 
in section 2.3, for the CME parameter /X 5 (t) in section 2.4, and, hnally, the magnetic field 
B{t) in the separate section 3. In section 4 we discuss our results and comment on some 
remaining problems. The detailed derivation of the induced electric current in eq. (1.1) is 
given in appendix A. 

2 Magnetic helicity evolution in a cooling NS core 

In this section we shall derive the system of kinetic equations for the continuous spectra of 
the helicity density and the magnetic energy density as well as for the chiral imbalance in 
the presence of the electroweak term V 5 in the current in eq. (1.1). We rewrite this system 
using the dimensionless variables in the form convenient for the numerical simulation, and 
obtain its solutions which obey certain initial conditions. 

We shall start with the Faraday equation for the magnetic field B evolution, 

— = a(t)VxB + 7?(t)V2B, ( 2 . 1 ) 

which can be obtained in a standard way from the Maxwell equations assuming the MHD 
approximation; cf. in refs. [12, 13]. Here a{t) = n(f)/crcond is the magnetic helicity parameter 
given by the anomalous current in eq. ( 1 . 1 ) with the coefficient 

20 - 

n(t) = ^(^ 5 (t) + F 5 ), ( 2 . 2 ) 

TT 

and r]{t) = (cTcond)”^ is the magnetic diffusion coefficient given by the electric conductivity 
o'cond- Note that one has to take into account the time dependence of cicond = c’'cond(i) in the 
cooling NS core; cf. ref. [14]. In deriving eq. (2.1), we suppose that the anomalous current 
in eq. (1.1) is added to the standard ohmic current Johm = o'cond(E + v x B) 

Basing on eq. (2.1), one can write the evolution equations for the binary combinations, 
h{t) ~ AB and psit) ~ B^. Here h{t) is the magnetic helicity density, 

h{t) = ^j d^x{A ■^) = yl • BD = / h{k,t)dk, (2.3) 

where V is the normalization volume, h{k,t) is the spectrum of the helicity density, and 
Bfc are the Fourier components of the vector potential A and the magnetic field B. The 
equation for h{k, t) should be complemented by the equation for the magnetic energy density 

^ (2.4) 

^Throughout the text we have neglected the bulk velocity evolution described by the Navier-Stokes equation 
since the length scale of the velocity variation is much shorter than the correlation distance of the magnetic 
field, A„ <C k~^. In other words, infrared modes of the magnetic field are practically unaffected by the velocity 
of plasma. In addition, the bulk velocity v does not contribute to the helicity evolution dh/dt ~ (E ■ B) when 
the generalized Ohm law is substituted, E = — v x B + 77V x B — aB. The small scale A„ is also a reason 
why we omitted the dynamo term V x (v x B) in the Faraday equation. Such a term could be important for 
a turbulent (early) stage of the NS evolution accounting for its differential rotation. 
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Note that, in eqs. (2.3) and (2.4), while dehning h{k,t) and p-Q{k,t) 


= ^^A{k,t) ■B*{k,t), pB{k,t) = ^^B{k,t)-B^k^t), (2.5) 

we perform the integration over the angles in the Fourier space meaning isotropic spectra. 

Accounting for the definition of h{k,t) and pB{k,t) in eq. (2.5), their evolution reads 
dth{k,t) ~ (Afc • + Afc • B^) and dtPBik,t) ~ (B^ • B^ + B^ • B^). Finally one gets that 

(see also Appendix D in ref. [15]): 


dh{k, t) 
dt ' 

dpB{k,t) 

dt 


2k‘^ 

^cond 

2k‘^ 


h{k,t) + 


(—) 

V '^cond / 




/ n \ 


<^cond 


PB{k,t) + - k^h{k,t). 

V <7cond / 


( 2 . 6 ) 


where 11 = n(t) is the parameter in eq. (2.2) responsible for the magnetic field instability. 
Note that analogous evolution equations were obtained in ref. [16]. 

Then we should derive the kinetic equation which governs the chiral imbalance P 5 {t) be¬ 
ing complementary to the system in eq. (2.6). For that purpose we use the Adler anomaly [17], 
— j'l) = i9/i(V’7^7^'i/^) = (2aem/7r)(E • B), and the magnetic helicity density evolution 


dh{t) 

dt 


2 

V 


/ 


d\(E-B), 


(2.7) 


which results from the Maxwell equations and eq. (2.3). Integrating the Adler anomaly, 
y f d^x(...), combined with eq (2.7), one gets the conservation law 


d 

dt 


nR - riL + 




TT 


-h(t) 


= 0 , 


( 2 . 8 ) 


where ur^l are the electron densities for right and left handed electrons given by their chemical 
potentials pr^l- 

Note that at the beginning of a chiral imbalance, pR ^ pL ^ p- Here p is the chemical 
potential of the ultrarelativistic degenerate electron gas which is fixed for the conventional 
abundance of electrons Ye ~ 0.05 in a cooling NS. For the nucleon density close to the nuclear 
one, riR ~ u-n = 0.18fm“^, the abundance Ye = UejnR = 0.05 corresponds to the electron 
density Ue = p^jSTr'^ = 9 x 10^®cm“^, or p = 125 MeV = const. 

Using the expression for the derivative of the density difference d(nR — nL)/dt ~ 
2p^{t)p^/ tt'^ and the hrst line in eq. (2.6) for dh{t)/dt = f dk[dh(k,t)/dt], the conserva¬ 
tion law in eq. (2.8) can be rewritten as the kinetic equation for the imbalance P 5 {t), 


dp^jt) 

dt 


/^^'^cond 


J dkk‘^h{k,t) 


dalmPBit) 


(Pbit) + Us) - Tfp^, 


(2.9) 


where we added the rate of chirality-flipping processes, Tf ^ {me/p)‘^i'co\h given by the 
Rutherford electron-proton (ep) collision frequency recoil = ^p/^^cond without flip. Here cOp = 
p^Aaem/^T^ is the plasma frequency in a degenerate ultrarelativistic electron gas. Note that 
the value of CTcond should correspond to a degenerate electron-proton plasma consisting of 
ultrarelativistic degenerate electrons and nonrelativistic degenerate protons. 
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Let us introduce the notations for the dimensionless functions: 


G? O? OL 

n{K,T) = n{K,T) = Mir) = (2.10) 

Z/J, /vminM 

where fcmin is the lower bound for the wave number k range (see below). Then, using eqs. (2.6) 
and (2.9), we get 


^ [-K^7i{K, t) + 2{M{t) + V)TI{k, t)] , 

^ [-k^TZ{k,t) + 2{M{t) + V)k^'H(k,t)] , 


dA^(r) 

dr 



K‘^ 7 i{K, T)dK 


2(7W(t) +V) 


TI{k, T)dK 


QM{t) 


.( 2 . 11 ) 


Here the argument k = fc/femin runs in the wave number region kmm = < k < femax, t = 

2 ^minV'^o is the dimensionless diffusion time, G = {aoTf/2k'^^^)/F = (2aem/3vr)(me/A:jnin)^ 
is the dimensionless rate of the chirality flip, and ao is the electric conductivity at the initial 
time to when the temperature of a NS core was Tq = 10® K. The factor F{t) = (To/ucond(t) = 
[r(r)/To]^ characterizes an increase of the electric conductivity (Tcond(t) ~ during the 
cooling of a NS core due to the neutrino (antineutrino) emission (see below). Finally, V = 
Q^emhs/'^^min = 7 X 10® is the dimensionless eN weak interaction potential for the fixed 
V 5 = 6 eV in eq. (1.2) given by the constant neutron density, Un 3> np, where Un = 0.18 fm“®. 

Note that the system in eq. (2.11) is a generalization of the kinetic equations derived 
in ref. [9] for the case of an arbitrary helicity density^ based on the system of the two kinetic 
eqs. (2.6) completed by the evolution of the chiral imbalance // 5 (f) in eq. (2.9). Indeed, 
considering the particular case of the monocromatic helicity density spectrum, T-L{k,t) = 
'H{t)6{k — kq), where 1 < kq < Kmaxj as well as assuming the constant conductivity with 
F = 1 and the maximal helicity density, at which 7^(r) = TI{t)/ko or h{t) = 2pB{t)/ko, we 
can recover the master equations in ref. [9]. 

There is, however, a discrepancy between eq. (2.11) and eq. (16) in ref. [9]. It consists 
in the additional factor 2 in the {Ai + V) term in eq. (2.11). This factor appears since in 
ref. [9] we have relied on the incorrect eq. ( 6 ) in ref [18]. The evolution of h{k,t) is correctly 
described by eq. (2.6) (see also refs. [15, 16]). Nevertheless, as we will see in sections 2.3-3, 
the main features of the magnetic field evolution, described in ref. [9], remain unchanged. 

Separating the magnetic diffusion factor ^diff('r) = exp F{T')dT''^ in the first 

two lines in eq. ( 2 . 11 ), one can easily find the important relation between the magnetic energy 
density spectrum TZ{k,t) and the magnetic helicity density spectrum 'H(k, r). 


^{k, t) 


[ 

TZ{k, t) 

2 

’77(K,ro)' 

/ 

K 


K 


(l-g2), 0 <<?<!. 


( 2 . 12 ) 


The relation in eq. (2.12) results from the conservation law for the auxiliary functions 
/A^mir) and TZi{k,t) = K^TZ{K,T)/Adis{T), 

^ [K^nfin, t) - njiK, r)] = 0, (2.13) 

which obey the simplified differential equations, drTiiinjT) = 2 F{t){M.{t) + V)TZi{k,t) and 
9 t-77i(k, r) = 2k^F{t){M.{t) +V)'Hi{k^t), obtained from eq. (2.11). 


- 5 - 

















At the initial time tq the relation in eq. (2.12) takes the form: 


n{K,To) = q 


K 


(2.14) 


Here the parametrization by the factor q < 1 corresponds to the the relation h{k,tQ) = 
2qp'B{k, to)/k in dimensional notations. Thns only the particular case q = 1 gives the maximal 
helicity density used in ref. [9]. 


2.1 Initial conditions 

In this section we choose the initial conditions for the system in eq. (2.11) corresponding to 
a realistic NS. 

The value of TZ{k,tq) = (aem/^min/.i^)/OB(fe) ^o) is given by the continuous initial mag¬ 
netic energy density spectrum [12, 13], 

PB{k,to) = Ck'^^"'^, femin < A: < /Cmax, (2.15) 


where the minimal wave number corresponds to the largest scale Ab = I?ns for an internal 
magnetic field within the NS core with the radius i?NS = 10 km, fcmin = = 2 X 10"^^ eV. 

The maximal wave number fcmax is a free parameter corresponding to the minimal spatial 
scale for the magnetic field A^™^ = /Cmax- The normalization constant C in eq. (2.15) results 
from eq. (2.4) and equals to 


(3 -|- v-b)Bq 

^^max 


(2.16) 


for fcmin ^ k^ascK- 

In eq. (2.15), we choose the Kolmogorov’s spectrum for the initial energy density with 
v-Q = —5/3, while other models of continuous spectra are possible, e.g. the Kazantsev’s 
spectrum with v-q = —1/2, or the white noise case i^b = 0. 

Basing on eqs. (2.15) and (2.16) and using eq. (2.14), we choose finally the following 
initial conditions in eq. (2.11) for dimensionless quantities: 


7i{K,To) = qC 


7^(K,To) = C 






1 < K < Kmax = 


kri 


» 1, 


C = 


Tmin 

«L(3 + z^b)So 


2p^k 


2 

max 


-Ad (to) = 


OLc, 


-pM = 1.2 X 10^^ 


Trkr, 


(2.17) 


Here in the last line P 5 {to) = IMeV <C p was used as an arbitrary value at the beginning 
of the chiral imbalance production through an electroweak mechanism at the earlier stages 
t < to of the NS evolution. 


2.2 Growfth of electric conductivity for a cooling non-superfiuid NS 

In this section we study the influence of the NS cooling on the time dependence of the 
conductivity and thus on the growth of the magnetic field. 

We remind that, in our problem, cicond is the electric conductivity in a degenerate 
electron-proton plasma consisting of ultrarelativistic degenerate electrons and non-relativistic 
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degenerate protons. The effects of both ee-collisions and the scattering of electrons by a neu¬ 
tron magnetic moment are minor for the electric conductivity in plasma, cicond = 
given mostly by ep-collisions [14], Here ujp = 4aem/37r is the plasma frequency in a degen¬ 
erate ultrarelativistic electron gas with the density rig = Note that in a degenerate 

electron gas both VcoW and Ucond depend on the temperature T since t'coii . This is due 

to the Pauli principle when all electron states with the electron momenta 0 < pe < h are 
busy, i.e. the ep-scattering is impossible at T = 0. 

For the constant electron abundance = ne/ne ~ const the electric conductivity rises 
with the cooling of a NS core [14]^, 

2 


'^cond(0 — *^0 


T{to) 

T{t) 




ao = 10 ' 


Up 


1036 


cm 


-3 


3/2 


MeV Ri const. 


(2.18) 


due to the negative derivative dT/dt < 0 given by the neutrino (antineutrino) emission [19, 

20 ]: 


dT(t) _ T{t) 
dt (riT — 2)t ’ 


(2.19) 


where the index ut equals to 6 or 8 depending on the neutrino emission channel'^. 

At early times, direct Urea processes, n ^ p + e~ + i^e and p -|- e“ —>■ n -|- t'e; prevail 
with the neutrino emissivity Ldu 10^^(r/10®K)®erg • s Thus, at this stage the index 
riT = 6 should be substituted in eq. (2.19). These equilibrium processes start just after a 
SN burst at first seconds, when the electron abundance is large enough Yp ~ 0.4. However, 
then the electron abundance Yp becomes less due to the deleptonization within the NS core 
as the SN remnant, and below Yp < 1/9 the direct Urea processes turn out to be suppressed 
because degenerate fermions are unable to conserve momentum while remaining on their 
Fermi surfaces. However, the presence of a third body, a nucleon N, in the modihed Urea 
processes, n-|-A^ —?> p-|-e“-|-A^-|-i>e andp-|-A^-|-e“ —)■ n-|-iV-|-t'e, allows to conserve momentum, 
and the neutrino (antineutrino) emissivity LmU ~ 10^'^(r/10® K)®erg • s“^, corresponding to 
the index nx = 8, dominates at Yp < 1/9. Thus we will adopt here nx = 8 in eq. (2.19) 
assuming also the conventional constant electron (proton) abundance Yp ~ 0.05. 

From eqs. (2.18) and (2.19) one can find the factor F = iTo/(Tcond(i) in eq. (2.11), 


' T{t) ' 
Tito)_ 


/ “2/(’^t-2) 

Vw 


( 2 . 20 ) 


or F{t) = (t/tq)”^^^ for nx = 8. Here we consider the cooling problem for a thermally 
relaxed non-superfluid NS core at the neutrino cooling stage. Before this time interval, during 
the first stage which lasts from ~ 10 yr to a few centuries, a newly born star is thermally 
non-relaxed. We put here to = 100 yr (Tq = 10® K) as the initial time (temperature) for 
our evolution eq. (2.11) and consider here the second stage which lasts for (10^ — 10®) yr, 

^The dependence of (Tcond in eq. (2.18) on the electron abundance, U = rie/riB, where ub = Un -t Up « 
is more slower than that on the temperature. For a cooling NS at the second stage the conventional value 
Yp ~ 0.05 corresponds to = 9 x 10^® cm“® 

^For the redshifted spatially constant internal temperature T{t) = r(r, t) exp[il?(r)], where T{r,t) is the 
local internal temperature, r is the radial coordinate and <E>(r) is the metric function that determines the 
gravitational redshift, we do not take into account in eq. (2.19) the constant factor exp[<l>(r)] « exp[<l>(i?Ns)] 
and consider the uniform internal temperature T{r,t) « Tb where Tb is temperature at the bottom of a thin 
envelope of NS related to the observable surface temperature by Ts oc \/Tb\ cf. ref. [20]. 
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to < t < 10® yr, when the core is thermally decoupled from the crust and it cools via the 
neutrino emission. During the third (final) stage, t > 10® yr, a thermally relaxed NS cools 
via the surface emission of thermal photons. Thus, the time to = 100 yr < t < 10® yr is 
appropriate in our numerical calculations where we use the cooling law in eq. (2.19). 


2.3 Magnetic helicity growth driven by weak eN interaction 


In this section we study the evolution of the magnetic helicity based on the numerical solution 
of the system in eq. (2.11). 

In figure 1 we show the corresponding growth of the helicity density h(t) = dkh{k, t) 
in the case of the Kolmogorov’s spectrum with = —5/3. Figures 1(a) and 1(b) correspond 
to fcmax = 2 X 10“^® eV or Ab = 1km, and figures 1(c) and 1(d) to fc max = 2 x 10~® eV or 
^(mm) _ ^ i?NS- To Support the helicity growth, i.e. to have dh/dt > 0, the maximal 

wave number fcmaxj that is equivalent to the smallest scale Ab = (fc m ax)~^, should obey the 
inequality®; 

kq < n(fo) = ^(/^ 5 (to) + ^ 5 ). (2.21) 

Z TT 

One can easily see the condition in eq. (2.21) is fulfilled for any chiral imbalance since 
Vs = 6eV = const S> kq for any k < /cmax- 

The most interesting issue here is the growth of the helicity density h{t) for the initially 
non-helical field, q = 0. It is remarkable that the difference of the helicity densities for the 
initially non-helical magnetic field, q = 0, shown by the solid lines, and the maximum helical 
field, q = 1 , shown by the dashed lines, vanishes at t S> to- In other words, an initial non- 
helical field tends to be maximum helical irrespective to the initial condition in eq. (2.17); cf. 
figures 1(a) and 1(c). On the contrary, at earlier times t >to such a difference is significant; 
cf. figures 1(b) and 1(d). In all cases we start from to = 100yr corresponding to the initial 
moment of the thermally relaxed non-superfluid NS core (see the comments in section 2.2). 

While the minimal wave number /cmin = is fixed, we can not choose rather big values 
of fcmax or the ratio Kmax = ^max/^min S> 1. Otherwise, the minimal scale of magnetic field 
Ajmm) _ could be comparable with a small scale of the fluid velocity we neglected 
here, when we consider large-scale magnetic fields and put A„ <C k~^, in order to avoid the 
involvement of the Navier-Stokes equation in our simplihed model. 

Considering smaller scales Ab = k~^, compared to those shown in lower panels in 
figure 1, we should also discuss earlier initial times to, when the stage of a thermally relaxed 
non-superfluid NS core have not been started yet. In other words, one expects that for a big 
^max the medium should be turbulent. This fact is not surprising because for small scales, 
i.e. when k is big, the evolution of h{t) and pB{t) proceeds faster since both characteristics, 
like the helicity density and the magnetic energy density, are proportional to the running 
wave number k, h{k,t) ~ Ak{t)Bk{t) ~ kA^ and pB{k,t) ~ involved in continuous 

spectra. The example of such an accelerated growth can be seen in figures 1(c) and 1(d), 
where fcmax = 2 x 10“® eV, or A^™^ = 100 m <C i?NS) is chosen. 

®Using eq. (2.6) one can show that the helicity growth is possible just at the initial time to if 


dh{k, t) 
dt 


t=to 


2pB(fc,fo) [_2gfe_^n(to)] >0. 

O-Q 
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Figure 1. The dimensional (in standard units G^cm) helicity density h{t) = 
{2fi'^kmax/ctlm) dK'H( k, t) versus time t' = t — to for the initially non-helical field, q = 0 

(solid lines), and the maximum helical field, q = 1 (dashed lines). Here fc„iin = = 2 x 10“^^ eV 

and to = 100yr. Panel (a) corresponds to fc^ax = 2 x 10“^° eV, or = 1km. The helicity 

evolution is shown for to < ^ < 4 x 10® yr. Panel (b) is the same as the panel (a) but for earlier times 
to < t < 5 X 10® yr. Panel (c) corresponds to fcmax = 2 x 10“®eV, or = 100 m ^ i?NS: for 

to < t < 10^ yr. Panel (d) is the same as the panel (c) but for earlier times to < t < 6 x 10^ yr. 

2.4 Change of the sign for CME 

In this section we discuss some peculiarities in the evolution of the chiral imbalance 

In figure 2 we show the evolution of the CME parameter A4 ~ /rs at different time 
scales based on the numerical solution of the system in eq. (2.11). In section 2.1 we have 
already commented on an initial positive ^ 5 (^ 0 ) ~ O(MeV) ^ /i = 125 MeV that vanishes 
fast during t ~ 10“®^ s <C to due to the huge chirality flip rate Q ~ 10®*^ [9]. This fast 
relaxation of /is to zero is not shown in figure 2. The example of the attenuation of /is —)• 0 
is given in ref. [9]. The curves in figure 2 start from the plateau A4 ~ 0. In hgure 2(b) one 
can see that this plateau lasts for ~ 2000 yr ;§> to for the maximal = 2 x 10~® eV. Then 
the CME parameter /is(t) changes sign, and for the same femax reaches the absolute value of 
the weak interaction parameter Vsj \-M.\ —)■ V = 7 x 10®, or |/is| —)• V5 = 6eV. Nevertheless, 
as seen in hgure 2, the sum Ai -\-V remains positive, A4 -|- V > 0, that provides the magnetic 
helicity growth, dTi/dr > 0, and the magnetic held strength itself, dTZ/dr > 0, as it results 
from eq. (2.11). 
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Figure 2. The dimensionless chiral imbalance A4 versus time t' = t — to for the initially non-helical 
held, q = 0 (solid line), and the maximum helical held, q = 1 (dashed line). We choose to = 100yr 
and fcmin = .^Ns = 2 X 10“^^ eV. Panel (a) corresponds to k^a.x = 2 x 10“^° eV, or = 1 km, and 

to < t < 4 X 10® yr. Panel (b) is built for = 2 x 10“®eV, or = 100 m, and to < t < 10^ yr. 


The change of the sign for the CME parameter /is at a certain moment of time, /is > 
0 —>• /is < 0 , is explained by the huge negative back reaction from the growing magnetic field, 
when the second term in eq. (2.9), ~ corresponding to the second negative term in the 

last line in eq. ( 2 . 11 ) for dA4/dr, becomes much greater than the chirality flip term, ~ Tj. 
For the initial magnetic field, Bq = 10^^ G, or at earlier times t < 2000 yr, the situation is 
reversal: the back reaction from a moderate magnetic field is negligible and the chirality flip 
term is the main one, which reduces the initial positive /is > 0 to zero, A4 —?• 0 . 

3 Generation of magnetic fields in magnetars driven by weak eN interac¬ 
tions 

The previous attempts to explain the growth of a seed field Bq ~ 10^^ G in a NS core up to 
the strongest B = (10^^ — 10^®) G observable in magnetars [4], which were based, e.g., on 
GME with /is 7 ^ 0 in refs. [21, 22], failed because of the underestimated chirality flip ~ Ey, 
see comments on this issue in refs. [9, 23-25]. 

In figure 3, we show that the magnetic held instability caused by the current in eq. (1.1) 
leads to the growth of a seed held Bq by about hve orders of magnitude. The resulting 
magnetic held grows up to H ~ 10^^ G for Bq = 10^^ G. The magnetic held growth happens 
in the interval ( 10 ^ — 10 ®) yr, depending on the minimal scale of f^e continuous 

spectrum in eq. (2.15). We revealed above the leading role of the weak interaction term V 5 for 
the magnetic held instability providing the growth of the magnetic held strength, dlZjdT > 0. 

Now it should be noted that our model (see also ref. [9]) does not require any special 
initial conditions, like an extremely strong magnetic held in a protostar in ref. [7], or a fast 
rotation in ref. [ 8 ], or a signihcant initial chiral imbalance in refs. [ 21 , 22 , 26], to generate 
strong magnetic helds in magnetars. Thus, using our results, we can predict a rather high 
magnetars abundance. This fact is in agreement with the hndings of ref. [27] that the 
number of magnetars in our universe should be comparable with that of ordinary pulsars 
having moderate magnetic helds Bq = 10^2 G 
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Figure 3. The dimensional (in Gauss) magnetic field in NS, B(t) = ^2 J dkpB(k,t), versus time 
t' = t — to for the initially non-helical field, q = 0 (solid lines), and the maximum helical field, 9 = 1 
(dashed lines). Here fcmin = = 2 x 10“^^ eV and to = 100yr. Analogously to figures 1 and 2, 

the Kolmogorov’s initial spectrum pB(k,to), given by eq. (2.15) with vb = —5/3, is assumed. Panel 
(a) corresponds to fcmax = 2 x 10“^° eV, or = 1km. The magnetic field evolution is shown for 

to < t < 4 X 10® yr. Panel (b) is the same as panel (a) but at earlier times to < t < 5 x 10® yr. Panel 

(c) corresponds to kjna.x = 2 x 10“®eV, or A|^'"^ = 100 m <C Rns, and times to < t < 10^ yr. Panel 

(d) is the same as the panel (c) but for earlier times to < t < 6 x 10^ yr. 

4 Discussion 

There are two anomalies in SM that give two corresponding terms in the electric current in 
eq. (1.1), J = JcME + Jcs- The former one is the Adler anomaly that leads to the well- 
known current in QCD plasma, Jcme = (2aem/7r)/r5B, see, e.g., ref. [28] or appendix A in 
ref. [29]. The latter anomaly leading to the Chern-Simons (CS) term Jcs ^ Gf arises due 
to the polarization effect when accounting for the parity violation in SM; cf. ref. [30]. The 
corresponding CS term in the SM Lagrangian Lcs = n 2 (A-B) = A-Jcs is given by the static 
limit of the parity violation term in the photon polarization operator Iljj ~ eijnknn 2 i^ = 
0, |k|) [31]. It arises as the second term in eq. (2.2) for large scale magnetic fields, |k| —)• 0, 
n 2 ( 0 , 0 ) = (2a 

em /ttM Gp. In the present work we did not apply the Matsubara technique 
to derive such 112, as it was made in our previous works [15, 32] for the electroweak ue- 
and ee-interactions. Instead we adopted a more transparent method to derive Jcs = II 2 B 
for massless electrons basing on the exact solution of the Dirac equation, accounting for the 
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electroweak eN interaction under the influence of an external magnetic field (see ref. [9] and 
Appendix A). This approach includes automatically the CME term Jcme in the current. 
Thus both terms, ~ and ~ V 5 , are present in the kinetic equations in eq. (2.11). The 
anomalous current J, being additive to the ohmic current enters Maxwell equation, results in 
the nonzero magnetic helicity parameter a = n/ucondj which governs the Faraday eq. (2.1) 
modified in SM due to anomalies. 

We generalized the approach for the magnetic field generation in a magnetar, which 
has been recently suggested in ref. [9], where the concept of the maximum helicity density 
was used. In the present work we considered an arbitrary initial magnetic helicity density 
in eq. (2.14), which is parametrized by the parameter 0 < g < 1. Here we also chose 
the continuous magnetic energy spectrum in eq. (2.15) instead of the monochromatic one in 
ref. [9]: pB{k, t) = pB{t)6{k — ko). The growth of the magnetic helicity density for the initially 
nonhelical field, h{k,to) = 0 corresponding to g = 0 , is owing to the presence of the term 
~ PB{k,t), which is nonzero at f = fo; in the first line in eq. (2.6). Although at early times 
nonhelical {q = 0 ) and maximum helical {q = 1 ) magnetic helicities are well distinguishable, 
the remarkable issue here is the tendency of the zero initial magnetic helicity [q = 0 ) to 
grow up to the maximal one {q = 1 ) at large time scales, compare figures 1 (a) and 1 (c) with 
figures 1(b) and 1(d). Of course, this helicity enhancement provides the similar dependence 
on the parameter q for the magnetic field itself in figure 3. 

As shown in section 2 . 2 , the cooling of a NS core due to the neutrino emission via the 
slow (modified) Urea processes leads to the increase of the electric conductivity. Such a 
cooling displaces with time the evolution curves for both h{t) and B{t) keeping the general 
effect of their growth at large times ( 10 ^ - 10 ^) yr. 

Note that our approach that is based on the consideration of the binary combinations, 
such as pb and h AB, gives no information on the structure and the orientation of the 
magnetic field. We deal here with random magnetic fields having Kolmogorov’s spectrum 
of the magnetic energy density. Such magnetic fields are characterized by the changing 
amplitude B and the varying large spatial scale Ab = k~^ in the continuous spectrum. 
Decreasing the minimal scale a|^™^ = k^l^^ down to the inhomogeneity size of velocities 
we should also consider the Navier-Stokes equation for the random fluid velocity v as well as 
distinguish a large scale magnetic field B and a small scale fluctuation b. Then the interplay 
of the kinetic helicity akin ~ (v- (V x v)), where the mean value is taken at large scales®, and 
the anomalous magnetic helicity parameter a = H/cTcond entering the Faraday eq. (2.1), both 
governing the evolution of large scale magnetic fields, would be important. In the present 
approach we do not take into account these problems. It should be also noted that in the case 
of the solid-state rotation with the spatially homogeneous angular velocity D(t) we assume 
here, the azimuthal rotation velocity of NS as a whole, Vq = D(t)rsin0e$, 0 < r < iiNSj 
does not contribute to the dynamo term V x (V X B), where V = Vq + v is the total fluid 
velocity and v is the small scale (random) velocity vanishing at large scales. As a result, the 
aU dynamo is not essential for us, while a^-dynamo is active, see, e.g., refs. [13, 15]. 

To resume we have further developed a novel mechanism, initially proposed in ref. [9], 
for the generation of the strongest magnetic helds in magnetars taking into account: (i) the 
neutrino cooling of the NS core, (ii) continuous magnetic helicity density and magnetic energy 
density spectra, and assuming (iii) an arbitrary initial magnetic helicity, including the case 
of the initially non-helical magnetic field. We have found that the CS anomaly caused by 

®In the dynamo term V x (v x b), the mean vector (v x b) « dkinB defines the pseudoscalar Okin ~ 
(v.(Vxv)). 
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the parity violation in the SM interaction of electrons with nucleons drives the growth of a 
seed magnetic field by about five orders of magnitude and its helicity by about ten orders 
of magnitude independently of an initial magnetic helicity density h{to) that could be even 
zero, h{to) = 0. Such a sharp magnetic helicity growth, followed by the strong enhancement 
of the magnetic field, is the main result of the present work. 
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A Averaged electron current induced by chiral effects 


In this appendix we obtain the solution of the Dirac equation for a massless electron inter¬ 
acting with nucleons under the influence of the magnetic field. Then we derive the induced 
electric current along the magnetic field direction. 

The Dirac equation for massless electrons reads 

[7^ (id^ + eA^) - 7° (UlPl + VnPn)] = 0 , (A.l) 

where 7 ^ = ( 7 ^, 7 ) are the Dirac matrices, = (0,0,Rx,0) is the vector potential for the 
magnetic held directed along the z-axis, Fl,r = ( 1 T 7^)/2 are the chiral projection operators, 
76 _ e > 0 is the absolute value of the electron charge, and Ul,r are the effective 

potentials for the interaction of left and right electrons with nucleons. The explicit form of 
hL,R can be found in ref. [9]. The upper sign both in =p and ±, throughout this appendix, 
corresponds to the left particles. 

Let us decompose V’e in the chiral projections as ipe = '4’h + V'Ri where V'l.r = -Pl.rV’e- 
Using the Dirac matrices in the standard representation [33], 


^ ' 0-1 


7 = 


0 cr 

—a 0 



(A.2) 


where cr are the Pauli matrices, it is convenient to represent '?/’lr = (v’l.Rj TV^l,r)- Using 
eq. (A.l) and separating the variables (/7 l,r = following 

equation for <^L,R(a^): 


[Po ± (crP)] V^L,R 


f Po±Pz TiVeP [d^ + ^] _ n 

V Ti\/^ [d^ -rj] PoTPz ) 


(A.3) 


where P^ = (Pl.r - 14,R, -^dx,Py + eBx,Pz) and rj = y/eBx + Py/y/eB. 

The solution of eq. (A.3) can be found using the Hermite function Un(??) = {eB 
X exp(—r 7 ^/ 2 )Ffn(r/)/\/ 2 '^n!, where is the Hermite polynomial, as 


¥^L,R(a;) 


1 




y/Po PzUn—l \ 
=FiV-Po ± PzUn ) 


(A.4) 
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for n = 1,2,.... The normalization coefficient in eq. (A.4) corresponds to the following 
normalization of the four component wave function V’lr ~ TV^l,r): 

/ (V'L,r)1p^p^ (V'L,R)„,p, p, {py - Py) S {p, - p',) . (A.5) 

The energy levels can be found from the expression, 

^0 = (^l,r - ^l,r)" =pI + 2 eBn. (A. 6 ) 


To obtain the solution in eq. (A.4) we use the following properties of the Hermite functions: 
[(9p + 77 ] Un = v^Un-i and [9p — 77 ] u^-i = — V^Un. 

If n = 0, the solution of eq. (A.3) has the form, 


<y5L,R(a:) 


1 

27rV2 



(A.7) 


Substituting eq. (A.7) to eq. (A.3) and using eq. (A. 6 ), we get that Pz > 0 for left particles 
and Pz < 0 for the right ones. It should be noted that, at n > 0, —00 < Pz < + 00 . 

Finally, we obtain the four component wave function in the form, 


/ V El,R - 14,R T PzUn-l \ 






U p 

V 

0 


1 


Ti\/ — 14, R ± PzUn 


T\/ A'l,R — Ul,R T PzUn-l 
\ i-\/ A'l,R — 14, R ± PzUn / 



( ° ^ 




1 

Uo 



2-Ky/2 

0 

\Two/ 




(A. 8 ) 


Using eq. (A. 8 ), we can calculate the averaged electric current along the magnetic field 


as 


r+oo 


/ dpy dpz [V’L,R7^V'L,R/e(^^L,R) “ V’L,R7VL,R/e(^^L,R)] > (A.9) 


n=0' 


where fe,e{E) = [exp{/3{E =p pl,r) + 1]~^ is the Fermi-Dirac distribution with upper (lower) 
sign ahead chemical potentials for electrons (positrons), = y^p^ + 2eBn ± Ul,r are the 
energy levels for electrons (upper sign) and positrons (lower sign), /3 = 1/T is the recipro¬ 
cal temperature, and pl,r is the chemical potential. The positron wave functions 7 /l£ in 
eq. (A.9) can be obtained from the electron ones t/^lR’ given in eq. (A. 8 ), by the charge 
conjugation [ 10 ]. 

At n > 0 we get for electrons 


'0L,R7^^L,R = T 


8 vr 2 (F;L R — Ul,r) 


X [(A'l,R - Ul,R T Pz) - (-EIl^ - Ul,R ± Pz) «n] 


and 


/ +00 

-00 


j 7 3 , eS p^ 

dpj/UL,R 7 Ul,r = -77—• 

Ul,r - Ul,r 


(A.IO) 

(A.ll) 
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Integrating eq. (A.ll) over one obtains that Landau levels with n > 0 do not contribute 
to the electric current. 

For the lowest Landau level with n = 0 we have 

r+oo _ 

''Al,r 7 ^V'l,r = J dpyV’L,R 7 ^V'L,R = (A. 12 ) 

Finally one obtains for left electrons 


g2 o ^+oo 

Jz = dpJeiPz + Ll), (A.13) 

and for right electrons 

j dpJe{-Pz + Lr), (A. 14) 

and the analogous expression for positrons. 

Adding the contribution of positrons with the properly changed signs, we obtain that 
the total current J = Jl + Jr induced by chiral effects reads 


Jz 



I [ <ipz [fe {-Pz + Lr) - fe {-Pz - 14)] 

L J —oo 

+ 00 'j 

[fe {Pz + Vl) - fe [Pz - Vl)] \ , 


or hnally (see eq. ( 1 . 1 ) above): 

J = ^(,.5 + U)B, 

TT 


(A.15) 


(A.16) 


which is additive to the ohmic current Johm in a standard QED plasma and where aem = 
e^/dvr is the fine structure constant, = (/tr — ^l)/ 2, and V 5 = (1/^ — 14)/2. It is worth 
to mention that eq. (A.16) is valid for any electron temperature. The first term ( Pb) 
in eq. (A.16) determines CME exploited, e.g., in QCD plasma [28], while the second term 
(~ V 5 ) is given by weak eN interactions in SM and has the polarization origin (compare in 
ref. [30]). 
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